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ABSTRACT 

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize 
performance and minimize costs. As a result, we have been investigating high-speed computational procedures to 
allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the 
progress towards a full three-dimensional and computationally efficient deterministic code for which the current 
HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann 
equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard 
Finite Element Method (FEM) geometry common to engineering design practice enabling development of 
integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 
milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential 
means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re- 
configurable computing and could be utilized in the final design as verification of the deterministic method 
optimized design. 


INTRODUCTION 

Following the Apollo mission, human exploration and development of space (HEDS) activity has been 
limited to low Earth orbit (LEO) with the advancing International Space Station (ISS) providing a first level of 
infrastructure for human exploration. Preparation for additional F1EDS activity beyond LEO is the natural extension 
of the current space program. Without dispute, providing protection for astronauts and equipment from the hazards 
of space radiation is a critical enabling technology for future HEDS activity and one of NASA’s two highest 
priorities (O’Keefe, 2002). Faced with a limited budget and an expanding space exploration program, the old way 
of doing business is inadequate and NASA requires revolutionary technologies to make advances. In none of the 
NASA enterprises is this more apparent than in HEDS. Radiation health risk mitigation will follow a triage 
involving biological, operational, and shielding countermeasures (Cucinotta et al. 200 1 ). Environmental modeling 
and shielding will play important roles in this development (Wilson et al. 1997-see especially Chapters 1 and 7, 
also see Wilson et al. 2001a). 

Early methods of space radiation shield evaluation relied largely on Monte Carlo codes (Alsmiller 1967, 
Lambiotte et al. 1971) and made important contributions to NASA engineering programs. Yet, slow computational 
procedures did not allow progress in coupling to spacecraft geometry for over 30 years (compare Alsmiller et al. 
1972, Annstrong et al. 2001) and did not allow early entry of radiation constraints into the design process and off- 
optimum Monte Carlo solutions to shielding problems continue to plague final designs (ASAP 2000, Qualls et al. 



2001, Wilson et al. 2001a). Simulations with full 3D Monte Carlo codes often use questionably simplified ID 
shielding geometry models to increase computational speed in which shielding for an ISS module is approximated 
as a very long aluminum cylindrical shell with 20.7 g/cm 2 thickness (Badhwar et al. 2001, Armstrong and Colbom 
2001) leading to an overestimate of the neutron flux and buildup of secondary charged particles within ISS since 
neutron leakage is underestimated (Clowdsley et al. 2002) and charged particle shielding is overestimated in these 
homogenized configurations (Wilson et al 1995a). Charged particle penetration estimates are improved in the ID 
Monte Carlo implementation of Pinsky et al. (2001) but neutron leakage is still incorrect. A detailed examination 
of the influence of these various ID simplified spacecraft geometries on the estimates of the interior environments 
are foimd to be large by Clowdsley et al. (2002). Furthermore the Monte Carlo implementations so far have been 
limited to a single construction material to maximize computational efficiency and optimization on multilayered 
optimized structures will not be possible using Monte Carlo methods because of greatly increased computational 
requirements. 

The “3D calculations” used in current ISS Monte Carlo studies are routinely exceeded in their correctness 
in coupling to spacecraft geometry by HZETRN (Wilson et al. 1995a, Shinn et al. 1998, Clowdsley et al. 2001a, 

2002, Qualls et al. 2001). In addition to spacecraft geometry simplifications, the astronauts in Monte Carlo 
calculations are approximated by a ID sphere of tissue. Conversely, the use of a 3D Monte Carlo code within such 
a simplified ID geometry for spacecraft and astronaut will often lead to erroneous underestimates of the solar 
particle event exposure since anisotropic shield distributions within an ISS like module have long been known to be 
a major factor (2 to 5) in specific astronaut organ exposure (Wilson et al. 1995a). The reason for such 
simplifications using Monte Carlo methods is easy to understand. To map out the interior environment of an ISS 
module requires at least 10° events within each volume about a field point. With 10’ points for a field map, this 
requires 10 9 events to map the interior fields. Each particle track from boundary to exit will experience 10 3 -10 4 
events requiring a ray trace per event leading to lO'MO 1 ’ ray traces for a field evaluation. In our ISS 16A 
configuration model in which only the F1AB module is currently represented in detail (Qualls et al. 2001), 14 
seconds is required to evaluate 10 ’ rays with a 12 processor script on a Onyx2 parallel machine or approximately 14 
msec per ray using optimized ray tracing methods (Qualls and Boykin 1997). The Monte Carlo interior field 
evaluation on a 12 processor Onyx2 would require lO '-IO seconds. It should be clear that Monte Carlo methods 
using ordinary central processor based computing technology will have limited usefulness in multidisciplinary 
design optimization. 

The development of high-speed deterministic procedures allows early entry of radiation constraints into the 
design process without simplification of the shield geometry (Wilson and Khandelwal 1974) but Monte Carlo 
methods can still play a role in future final design evaluation with full geometry provided more efficient ray trace 
methods can be found, massively parallel machines are used, and/or a radically different approach to ray tracing can 
be found. We propose a radically new approach based on the use of the emerging re-configurable computing 
technology to enable Monte Carlo methods for design evaluation (but not optimization) and at the same time allow 
deterministic methods to map the ISS interior fields in 10’s of minutes (near real-time evaluation). 

At this juncture, the deterministic evaluation of multiple-elastic scattering of heavy ions has been 
adequately resolved with detailed laboratory testing (Schimmerling et al 1986, Shavers et al. 1990, 1993). The 
perplexity of straggling has been resolved (Wilson and Tai 2000, Wilson et al. 2001b, 2002). The 3D geometry in 
charged particle transport can be accurately approximated by added asymptotic terms (Wilson and Khandelwal 
1974, Wilson 1977). The lowest order asymptotic tenn has been solved as a marching procedure (Wilson and 
Badavi 1986), a perturbation expansion (Wilson et al. 1984), a non-perturbative expansion (Wilson et al 1994a), 
with extensive testing in laboratory and flight experiments. Detailed investigations of the diffusive neutron 
component at various levels of approximation from 3D Sn methods (Singleterry and Wilson 1998) to collocation 
methods (Clowdsley et al. 2000) including PN methods have shown promise in approaching an efficient full 3D 
code but needs further development. It has been demonstrated that these methods can be used efficiently in the full 
3D complexity of spacecraft design and have been flight-tested (Badhwar et al. 1995, Shinn et al. 1998, Clowdsley 
et al. 2002). Added treatment of the charged particle diffusive components, which will converge more rapidly than 
the neutral neutrons can be adequately treated using perturbation theory. A meson database needs completion 
(Cucinotta et al. 1998, Norbury et al. 2002) with addition of a complete E&M cascade code (Nealy et al. 2002). 
Further considerations of fast computational procedures and implementation issues on high performance computers 
for ultimate insertion into shield optimization and reliability based shield design methods are now being made 
(Wilson et al. 2001a, 2002, Tripathi et al. 2001, 2002). 

In the present report, we will discuss the development of a fully 3D deterministic code for which past code 
development resulting in F1ZETRN is the lowest order approximation. We will also discuss the requirements for 



implementing a hypercomputer algorithm to enable a potentially high-performance Monte Carlo version of LAHET 
useful in spacecraft design evaluation. 

DETERMINISTIC CODE DEVELOPMENT 


The relevant transport equations are the linear Boltzmann equations derived on the basis of conservation 
principles (Wilson et al. 1991a) for the flux density (f>j(x, Q, E) for particle type j as 

£}• V<fij(x,£2,E)= TJ Ojk(Qfi’’E,E) <f>k(x,£2’,E) d£2’ dE' - Oj(E) (f)j(x,£2,E) (1) 

where <jj(E) and Ojk(Q,£2, E,E) are the shield media macroscopic cross sections. The <7j -±(£2,0:' ,E,E) represent all 
those processes by which type k particles moving in direction Q' with energy E' produce a type j particle in 
direction £2 with energy E ( including decay processes). Note that there may be several reactions that produce a 
particular product, and the appropriate cross sections for equation (1) are the inclusive ones. Exclusive processes 
are functions of the particle fields and may be included once the particle fields are known. The total cross section 
(7j (E) with the medium for each particle type is 


<5j (E) = (Jj >at (E) + (7j ie l (E) + Oj r (E) 


( 2 ) 


where the first tenn refers to collision with atomic electrons, the second term is for elastic nuclear scattering, and 
the third term describes nuclear reactions where we have ignored the minor nuclear inelastic processes. The 
corresponding differential cross section is similar. Many atomic collisions (~ 10 6 ) occur in a centimeter of ordinary 
matter, whereas ~ 10 3 nuclear coulomb elastic collisions occur per centimeter, while nuclear reactions are separated 
by a fraction to many centimeters depending on energy and particle type. Solution methods use the atomic 
collisions as a first order physical perturbation with special methods used for neutrons for which atomic cross- 
sections are zero. 

The first order physical perturbation to 
the right side of equation (1) is the 
atomic/molecular cross sections as noted in 
equation (2) for which those terms in (1) are 
expanded about the energy moments resulting in 
range/energy relations including relativistic 
polarization effects (Shinn et al. 1993) and 
straggling parameters (Wilson et al. 1991, 

Wilson and Tai 2000, Wilson et al. 2001b, 

2002). The distribution of the electrons about 
the ion path is critical to evaluation of biological 
injury (Cucinotta et al. 1995), is critical to the 
evaluation of shield attenuation properties 
(Wilson et al. 1995b), and fundamental to 
dosimetric evaluation of astronaut exposure 
risks (Shinn et al 1999). Such effects are 
likewise governed by equation (1). The next 
physical perturbation term is the coulomb 
scattering by the atomic nucleus and is 
represented by Rutherford scattering modified 
by screening of the nuclear charge by the orbital 
electrons using the Thomas-Fenni distribution 
for the atomic orbitals. The total nuclear coulomb cross section found by integrating over the scattering directions 
is related to the radiation length (Schimmerling et al. 1986) and has been well validated in laboratory tests (Shavers 
et al. 1990, 1993). The electron production in ion collisions and the electromagnetic cascade following the neutral 
meson decay are described, in part, elsewhere (Cucinotta et al. 1998, Anderson et al. 2001, Nealy et al. 2002). 

The nuclear reactive differential cross sections can be written in the following fonn 



Fig. 1 Isotropic and forward neutron spectra produced by 500 
MeV proton in aluminum. 


(j jk J£2,£2’,E,E ) = Oj kiso (E,E ) + cr jkfor (£2.£2\E,E ) (3) 

where the first tenn is isotropic and associated with lower energy particles produced including target fragments and 
the second tenn is highly peaked in the forward direction and is associated mainly with direct quasi-elastic events 



and projectile fragments (Wilson et al. 1988). Surprisingly even nucleon-induced reactions follow this simple form 
and the isotropic term extends to relatively high energies (see Fig. 1). For nucleon induced reactions, the following 
form has been used in versions of FLUKA (Ranft 1980) as follows 

Ofr/aa&E) = v jk (E) <j jkr (E )J jk (E,E ) g R (Q»Q,EA T ) (4) 

where the Ranft factor used in FLUKA is 

g R (£2»i2,E,A T ) = N R exp[-^/X R ] jt/2>0>0 ( 5 ) 

and constant for larger values of production angle #and X R given by Ranft as 

X R = (0.12 + 0.00036A t /E) (6) 

although new generalized fits are being derived. This separation in phase space will be further exploited in 
computational procedures. The heavy ion projectile fragment cross sections are further represented by 

a jKfor (Q,Q’,E,E) = Ojk./E) N, exp [-2m V(EE ) ( 1-Q»£i’)/ e t j k ] exp[- (E + /L jk -E'f/2 e jk J/V(2; Z£ jk ) (7) 


where X jk is related to the momentum downshift, e jk 
is related to the longitudinal momentum width, £ tJk is 
related to the transverse momentum width, and N, is 
the transverse normalizing factor. Since the 
transverse width is small compared to the projectile 
and fragment energy the transverse function is highly 
peaked about the forward direction (£2»Q' ~ 1). 

Atomic interactions limit the contributions of 
charged particles in the transport process. For 
example, the protons and alpha particles produced in 
aluminum below 100 A MeV contribute to the 
fluence only within a few centimeters of their 
collision source and the heavier ions are even more 
restricted (see Fig. 2). This is an important factor in 
that the transported secondary charged particle flux 
tends to be small at low energies and the role of 
additional nuclear reactions are likewise limited (see 
Fig. 3). 

We rewrite equation (1) in operator notation 


by defining a vector array field function as 

0 = [</>j(x, a E)] (8) 

the drift operator 

D = [n»Vj (9) 

and the interaction operator 

I /Tj <jj[(£2,Q', EE) d£2' clE’ - <jj(E)J (10) 

with the understanding that / has three parts associated 
with atomic, elastic, and reactive processes as given in 
equation (2). Equation (1) is then rewritten as 

[D-l at -l el ]»0 l r »0 (11) 


where the first two physical perturbation terms are 
shown on the left-hand side and have been adequately 
resolved for ions in past research. The reaction cross 
section is separated by equation (3) into isotropic and 
forward component for which equation (11) may be 
written as coupled equations 



Fig. 2 Range of ions in aluminum. 



Fig. 3 Probability of nuclear reaction as a function of ion 
type and energy. 



[D - I at -I el + <3,1 • 0for = {\ O rf JQ,Q,E,E') dQ dE'}» 0 for =S rfor & for (12) 

and 

[D-I.,-I d + <T r ]»0 iso = t\ a(Q,£i',E,E) da dE'j* 0 is „+ t\<7 r .jnj7,E,E) dil dE’!*0 for =E r . 0 isl> + S^„» (13) 

Equation (12) can be written as a Volterra equation (Wilson 1977, Wilson et al. 1991) and solved as 

0for =[G+ G»S, Jo ,»G + G»S ly/or »G»S r/or »G (14) 

for which the series can be either evaluated directly or proscribed as a marching procedure in either a perturbative 
sense as the current form of HZETRN or nonperturbative sense (future version of HZETRN) as described 
elsewhere (Wilson et al. 1994a). 

The cross term in equation (13) gives rise to an isotropic source of light ions and neutrons of only modest 
energies for which Fig. 1 is typical. The high-energy portion of the isotropic spectra arises from multiple scattering 
effects and the Fermi motion of the struck nucleons within the nucleus. The low-energy isotropic spectra arise from 
nuclear decay processes. Spectral contributions to the Neuman series depend on the particle range and probability 
of surviving nuclear reactions that establish the functional form of the G matrix. The second term of the Neuman 
series is proportional to the probability of nuclear reaction that is limited by the particle range as discussed above 
and shown in Fig. 3. It is clear from Fig. 3 that those nuclear reactions for the charged particles below a few 
hundred A MeV are infrequent for which fast convergence is expected. For the moment we will neglect the 
straggling and multiple-elastic processes to simplify the present explanation (provide only minor corrections to 
space radiation exposures but important in laboratory testing) and examine the remaining reactive terms of equation 
(14). The corresponding Volterra equation is given by 

<t>j(x,Q,E) = (Sj(E)Pj(E) §(H(npQ,Q£) + l/ E Er dE’A j P j (E)f E -LdE"d£2’ <j jkJ .JQQ,E’,E”) 

X </> k (x+[Rj(E) - R/E )jan’,E ■)}/ Sj(E)Pj(E) (15) 

where F is the point on the boundary connected to „v along -Q, E r = Rj 1 [p- d + RJ, p is the projection of x onto Q, 
and d is the projection of T onto Q. Equation (14) results from the Neuman series solution to equation (15). In the 
past we have expanded the angular integral Q asymptotically and implemented as a marching procedure 
(HZETRN, Wilson and Badavi 1986), as a perturbation expansion (Wilson et al. 1984), and by non-perturbative 
approximation (Wilson et al. 1994a) resulting in three distinct methods to evaluate the first order asymptotic terms, 
all of which have had extensive experimental validation (Shavers et al. 1993, Wilson et al. 1998, Shinn et al. 1998). 
Independent of the method used to evaluate the lowest order term, the first correction term is found by replacing the 
fluence in the integrand of equation (16) by the lowest order asymptotic solution as 

<t>j(x,n,E) = {Sj(E)Pj(EJ <i>j(r(Q r x),Q£) + lI E Er dE’A ] P ] (E)f E m LdE"dQ l (T jkf JQQ,E’,E”) 

X0Jx+[Rj(E) -Rj(E )]Q,Q£')}/ Sj(E)Pj(E) (16) 

where </>j(x,f2,l:) is found as an integral over the neighborhood of rays centered on Q using the lowest order 
asymptotic solution <f> ko (x,£2’,E”) along an adjacent ray 
directed along Q’. Note that the boundary condition reached 
along -Q’ enters through the lowest order asymptotic 
approximation and the angular integral correction in equation 
(16) is determined by the homogeneity and angular 
dependence of the space radiation and radius of curvature of 
the bounding material as we have shown long ago (Wilson and 
Khandelwal 1974). These are the determinant factors of the 
magnitude of the first order asymptotic correction which is 
anticipated to be very small for human rated systems (large 
radius of curvature) in space radiation which is homogeneous 
and isotropic in most applications (Wilson et al. 1991, Wilson 
et al. 1994b). 

In a region of small radius of curvature the specific 
flux components near the site of evaluation will be missing 
contributions along adjacent rays which do not compensate 
losses along the ray on which the solution is evaluated 
representing the losses due to leakage. (Note, an asymptotic 
treatment of such small angle dependent phenomena is the 
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Fig. 4 Normalized transverse components for 
Ca fragmentation. 




only useful approach circumventing large discretization errors.) This computational procedure is only a small 
addition to prior code development and will have little impact on computational efficiency. The angular 
dependence of the integral kernel of equation (16) is controlled by the forward reactive cross section 
<j jkfm .(£2,{2',E’,E”) with its highly peaked structure given by equations (4) or (7) depending on particle type. The 
angular dependence of the forward peak of fragmenting Ca ions at 100 and 1,000 A MeV is shown in Fig. 4. The 
low-energy ions with limited range have transverse components on the order of 10 degrees reducing to a few 
degrees at high energies. Note that the low energy ions have limited range and will contribute little to the 
transported flux (see Fig. 2) or nuclear reactions (see Fig. 3). The higher energy ions with their much longer 
pathlengths giving more important contributions are related to only a very small angle of acceptance (few degrees) 
at the boundary. The form of the kernel leads directly to a Gauss-Hermite expansion and evaluation over the angle 
of production. Although the neutron Neuman series for the forward components converge more slowly since their 
contribution to the neutron flux is not limited by atomic interactions these higher energy neutrons will be 
adequately evaluated by similar procedures. Fligher order asymptotic terms can be evaluated with similar iteration 
of equation (16) if required but all indications are that the first such correction will be small for space radiation. 
This leaves the diffuse components of neutrons and light ions produced in the collision of the forward components 
and transported by equation (13) to be resolved. 

The transport from the low-energy neutron and light -ion isotropic sources in equation (13) dominates the 
solution below about 70 A MeV (see Fig. 1). In this region light-ion transport is completely dominated by the 
atomic interaction terms and only a very small fraction have nuclear reactions making only minor contributions to 
the particle fields. This is especially true for the target fragments that can be solved in closed form (Wilson 1977, 
Cucinotta et al. 1991). The transport solution for the isotropic ion source terms to the lowest order perturbation is 
given by 

</>j >iso 0 (x,£2£) = lk Er dE ’A jP j(E )f E ™J 4k dE”d£2 ’ <T jkisor (£2,£I,E’,E”) </> kf Jx+ [R/E) - Rj(E’)]Qi2’,E ”)/ S/E)Pj(E ) 
(17) 

and will give highly accurate solutions to equation (13) since very few of the ions will have reactions (see Fig. 3) 
but could be easily corrected using the HZETRN light-ion propagator applied to the diffusive source terms. Note 
the E’ integral effectively sums the ion source terms along direction D from the boundary to x. Also the nuclear 
survival terms Pj(E) are all near unity (see Fig. 3 showing 1 - Pj(E)). 

The neutrons have no charge and are undergoing, at low-energies, mainly elastic but also reactive nuclear 
processes. Although, equation (13) exhibits behavior similar to thermal diffusion there are strong differences 
between thermal and neutron diffusive processes. Thermal diffusion at ordinary temperatures has minor leakage 
through near boundaries since radiative processes are proportional to T 4 (in the absence of convection) leaving 
lateral diffusion an important process. In distinction, neutron diffusion is dominated by leakage at near forward and 
backward boundaries and lateral diffusion plays a minor role. Generally, low-energy neutron leakage is a dominant 
process within 15-20 g/cm 2 of the bounding surface in most materials. Since human rated systems have shields of 
large radius of curvature and small thickness to radius ratio as determined by living and working space 
requirements, it approximates a connected system of flat plates for which leakage at forward and backward 
boundaries dominates the transport. In this limit, neutron transport simplifies to a sequence of approximate ID- 
transport problems with leakage at the back and forward boundaries and shows reasonable success in comparison 
with experimental flight data (Clowdsley et al. 2000, 2001). hi the present development we will consider a 
convergent series of approximations to gauge accuracy of the transport procedures and allow choices of the most 
practical method. 

The first term for diffusive neutron transport uses the lowest order perturbation similar to equation (17) 

given as 

<t> nMo °(x,OE) = Zf 0 p - d dx’exp[- o n (E)x] f E r/ 4x dE”d£2 ’ (T jkisor «2,a',E’,E”) 0 k/ Jx -x'QS2’,E’’J (18) 


where p-d is the distance along -Q from x to the boundary, x ’ is the distance from x to the source point along -£2, 
and cr„(E) is the total neutron cross section. Note that equations (17) and (18) transport all particles associated with 
the collisions of the forward component & for found as solution to equation (14). The second collision term 
associated with the diffuse charged particle field given by equation (17) is negligible but additional source terms 
from the lowest-order diffuse neutron solution given by equation (18) provides a strong source of diffuse neutrons. 
The added transport of the neutrons is given by 




where <j> njso 1 (x,QE) is the remaining diffuse neutron component. The source term to the far right of equation (19) 
have been solved in exact 3D geometry and the energy spectrum is much degraded for the source term in equation 
(18). Typical spectra of the sources in equations (18) and (19) are shown in Fig. 5. It is clear from the properties of 
the second collision source term that the diffuse spectra of neutrons from this term are highly degraded in energy 
and the methods developed in the nuclear engineering community for reactor applications are fully applicable and 
our final attention turns to solution of equation (19). Note that the software prepared for equation (19) is also 
applicable to nuclear reactor shielding issues with appropriate source terms. 

The dominant contribution to the low-energy neutron transport in most materials is elastic scattering from 
the media nuclei. In our shell FEM geometry, equation (19) can be written as 


m + <J n (E)]dJ(x,ii,E) = \ a<Q,£2\E,E)dQ\lE'd l J(x^',E') + \o n (QQ’,E,E’)dQ'dE'6 


(20) 


where the last term on the right is the source term. The elastic scattering term in equation (20) has a unique angular 
contribution jU 0 = £2»£2' for a given energy transfer (E’ - E) relating to a unique direction p’ under the integral. Standard 
spherical harmonic expansions of the cross sections and 
ilux are made to develop a solution. As example, first 
order Legendre expansion (PI approximation) gives the 
forward and backward llux as 


0f,b — 0-5 c^o - 0-75 01 


( 21 ) 
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where 0 O , 0i are the Legendre coefficients of <j> n t J(x,n,E) 
with approximate transport equation given as we have 
used in the past (Clowdsley et al. 2000) using multigroup 
and collocation methods. 

The multigroup transport of the diffuse 
neutrons in multilayered flat plate geometry with 
variable front and back boundaries in which the ions 
are treated in the lowest order asymptotic 
approximation is the current production HZETRN 
code. The first improvement would be to treat the 
diffuse ion components with perturbations from the 
neutron diffuse component giving a complete flat 
plate code in lowest asymptotic order with 
multigroup neutrons and diffuse ion components. To 
this we will add the first asymptotic correction. This 
will provide NASA with an interim HZETRN 
production code for engineering design process 
development. At each step of future improvements, 

the corresponding codes will be integrated into the SIREST/FACE collaborative engineering environment 
(Singleterry et al. 2001) with environmental models for NASA use (the current production code resides there 
already). To this code we will add mesons and electromagnetic cascades in the near future. 
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Fig. 5 Spectral dependence of diffuse isotropic neutron 
source term (eq. 18) and the first collision source term 
(eq. 19). 


MONTE CARLO/SPACECRAFT INTERFACE DEVELOPMENT 


Ray tracing procedures in combinatorial geometry is the traditional method of Monte Carlo and will not 
meet NASA requirements in that the models are difficult to build and don’t interface with other engineering 
analysis packages. The ray tracing procedures in current use with deterministic codes, built on CAD geometry 
methods, have their origin in methods for 2D-video display. To accomplish this, the CAD models are transferred to 
a faceted or FEM geometry representation and mapped into a WAVEFRONT format that allows rapid ray analysis 
(Qualls and Boykin 1997). A FORTRAN code XRADICAL was written to quickly analyze each facet and 
calculate intersections. XRADICAL procedures were designed to emulate on a serial machine (or with limited 
parallelization) the processing similar to that found using gate arrays on high-speed video cards. Standards on 
CAD model development to arrive at meaningful results have been given (Qualls and Boykin 1997). Although 
highly optimized, the limitation of the method is the calculation and testing of many facets in typical spacecraft 
applications in a serial or near serial manor. For example, the 14 msec per trace using 12 parallel processors on an 
ORIGIN machine in the simplified ISS geometry with 62,000 facets (only the HAB module has fine detail). 

The approach to resolve this computational task for Monte Carlo applications is to migrate back to gate 
array processing with large inherent parallel computational paths to improve the ray trace performance. To 



accomplish this task, we will investigate setting up ray trace procedures on a radically new computational platform 
(re-configurable computer), a field programmable gate array (FPGA) with a reasonably high level language 
developed by STARBRIDGE. This is a fast developing technology and Langley Research Center is fortunate to 
have an R&D machine and close working relationship with STARBRIDGE under a Space Act Agreement. The 
research under this element will be to develop ray tracing procedures for relatively simple objects but more 
complicated than current Monte Carlo can reasonably handle as a test of the use of this technology. The success 
will be determined on comparisons between the serial code and the FPGA supported code. The advantage of this 
computer is the nearly unbound degree of parallelization giving potentially high performance similar to the high 
performance video card implementation. Plans will be made to scale to more realistic geometries and alternate 
pathways to implement such a procedure on the programmable STARBRIDGE hypercomputer. In the later stages 
of this work, the interface between the hypercomputer and Monte Carlo codes on a serial machine will be 
considered. Preliminary tests using a serial machine will be developed and the whole system optimized. 

CONCLUSIONS 

The NASA approach to risk mitigation involves a triage of science -based methods and relies heavily on the 
development of an understanding of the initial physical insult and radiobiological response mechanisms. The first 
line of mitigation is control of the physical insult through shielding material arrangement in the domains of human 
occupancy. The approaches to this problem have been discussed in terms of current technical development and 
limitations. The pathways to a final resolution of a high-efficiency shield code for implementation into spacecraft 
design processes is clear requiring only development and validation. Indeed, a reasonable software and database 
has already been assembled and tested in both laboratory and space flight. Additional development will be in the 
main only minor improvements to the existing technology and validation. The exposure of human beings to 
personal health risks demands the highest level of ethical awareness in order to achieve exploration without 
compromising the respect due to the humanity of the explorers. 
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